A new strategy for probing the Majorana neutrino CP violating phases and masses 



David Delepine and Vannia Gonzalez Macias 
Departamento de Fisica, Universidad de Guanajuato, 
Campus Leon, CP. 37150, Leon, Guanajuato, Mexico. 

Shaaban Khalil 
Centre for Theoretical Physics, The British University in Egypt, 
El Sherouk City, Postal No, 11837, P.O. Box 43, Egypt 



o 
o 

(N 

< 



in: 



> 

oo 
in 

(N 
00 

O 

o\ 
o 



X 



Gabriel Lopez Castro 
Departamento de Fisica, Cinvestav, Apartado Postal 14-740, 07000 Mexico D.F., Mexico 

(Dated: August 15, 2009) 

We propose a new strategy for detecting the CP-violating phases and the effective mass of muon 
Majorana neutrinos by measuring observables associated with neutrino-antineutrino oscillations in 
TT* decays. Within the generic framework of quantum field theory, we compute the non-factorizable 
probability for producing a pair of same-charged muons in vr* decays as a distinctive signature of 
f/i — i^fi oscillations. We show that an intense neutrino beam through a long baseline experiment 
is favored for probing the Majorana phases. Using the neutrino-antineutrino oscillation probability 
reported by MINOS collaboration, a new stringent bound on the effective muon-neutrino mass is 
derived. 
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I. INTRODUCTION. 

Detecting CP violating phases in the lepton sector is 
one of the most challenging problems in the study of 
neutrino mixing. In the basis of diagonal charged lep- 
ton mass matrix, the neutrino mass matrix rui, can be 
written in flavor basis as — U*m'^^^^U\ where U 
is the PMNS mixing matrix, which can be written as 
U — Kdiag (l, e'"2, 6*"=*). Ui are the Majorana phases 
and the mixing matrix V can be parameterized by one 
Dirac phase and three angles: solar angle 0i2, atmo- 
spheric angle ^23, and Chooz angle 6113. 

The measurement of the effective electron-neutrino 
mass in the neutrinoless double beta decay (Oz^/3/3) exper- 
iments can not restrict the two Majorana CP violating 
phases present in the PMNS mixing matrix. The effective 
electron- neutrino mass {m^e) is given by 



\{mee)\ = \ Y.Ul 



(1) 



This effective mass parameter depends on the angles 
6*12 and 6113, the neutrino masses m^., Dirac CP phase, 
and Majorana phases a^. There are several studies on us- 
ing the results of (Oi^/3/3) together with the new data from 
terrestrial and astrophysical observation in order to re- 
strict the Majorana neutrino CP violating phases [l|, 
However, this analysis is model dependent and quite sen- 
sitive to the ansatz of the neutrino mass spectrum: quasi- 
degenerate, normal or inverted hierarchies. In this re- 
spect, it is not possible to measure the Majorana neu- 
trino CP phases from (Oi/(3(3) experiment. This may be 
expected since in the {Qv(30) one measures the lifetime of 



the decay of two neutrons in a nucleus into two protons 
and two electrons, which is a CP conserving quantity. 

On the other hand, direct bounds on other effec- 
tive neutrino mass parameters (to/;) = ^iUfimi^ from 
present experimental data are very poor. Currently, the 
strongest bound for the muon-neutrino case from the 
fi'^ 11'^ branching fraction is only {mf^f^) < 
0.04 TeV [5], which leads to a negligible constraint on 
the neutrino masses and CP violating phases. Therefore, 
it is commonly believed that direct bounds from other 
AL = 2 decays are only of academic interest and can not 
fix the neutrino mixing parameters Q . Some attempts to 
detect CP violation based on the difference between os- 
cillation probabilities of neutrinos and antineutrinos can 
be found in Ref.0]- 

Here we propose a mechanism, based on neutrino- 
antineutrino oscillation [1, [^, [l3| , which would allow to 
derive a strong bound on the effective mass of the muon- 
neutrino. In addition, it provides a method for detect- 
ing the Majorana neutrino CP violating phases through 
measuring the CP asymmetry of the tt^ decay where 
neutrino-antineutrino oscillation take place. Using the 
bound on the neutrino-antineutrino oscillation probabil- 
ity reported by the MINOS Collaboration [11], we derive 
a bound on (m^^) which improves existing bounds by 
several orders of magnitude. 

It is worth noting that the probability of a process 
associated to neutrino oscillation is usually assumed to 
be factorized into three independent parts: the produc- 
tion process, the oscillation probability and the detec- 



tion cross section. In Ref. I4I, this approximation was 



avoided and a generic framework based on quantum field 
theory was proposed to get a simple expression for the 
CP asymmetry. Here, we adopt the S-matrix amplitude 
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FIG. 1: Feynman diagram of the process it (pi) ^i^(p2) 
i^/j(p) followed by the detection process:7^(p) + N{pn) 
N'{p^.)+li+{pi). 



method described in lli 



II. NEUTRINO-ANTINEUTRINO 
OSCILLATION 

Let us start by considering a positive charged pion 
which decays into a virtual neutrino at the space-time 
location (x, t) together with a positive charged muon. 
After propagating, the neutrino can be converted into an 
antineutrino which produces a positive charged muon at 
the point {x ,t) when it interacts with a target, as shown 
in Fig. 1. For definiteness, we illustrate this process with 
the production of the neutrino in 7r+ decay and its later 
detection via its weak interaction with a target nucleon 
N 



^ V-i\p) + N{pn) 



N'{p 



N' , 



where the superscript s{d) refers to the virtual neu- 
trino (antineutrino) at the source (detection) vertex. 
This |AL| = 2 process is a clear signal for neutrino- 
antineutrino oscillations of the muon type and its am- 
plitude should be proportionnal to neutrino Majorana 
masses. 

If one ignores other flavors, the time evolution of the 
v^i — v^ system would be analogous to that of the K^ — K^ 
or — systems. Instead, we prefer to use the for- 
malism developed in Ref . . where the whole reaction 
includes the production and detection processes of neu- 
trinos. The decay amplitude becomes (for simplicity we 
assume that leptonic flavor is conserved at the production 
and detection vertices): 



^{G pVudY [J N N') l^fT^ 

X 5I^A'(p;)7^(i + 75) i)iv{p2) 



-HE,, 



(2) 



where the relation Vk = '^Ukai^a between flavor k and 
mass a neutrino eigenstates has been used, /tt = 130.4 



MeV is the tt^ decay constant, and J^^v' parametrizes 
the interaction with the nucleon. Note that, contrary to 
the case of neutrino oscillations ,12], only the neutrino 
mass term survives in this case. For simplicity, one as- 
sumes that 

{JNN')i^ ^Uj^'{Pf^')jf,{gv+gAl5)uN{PN) (3) 

where we use qy = gvil^ = 0) = 1 and gA = gAil^ — 
0) w —1.27 [ij]. If we neglect terms of 0(m^/mjY jy')j 
one obtains 



= i2n)^S\pi-pN+PN' +P2-pi)iGFV^dr 
X %^{gA-\fraNTal{E2-E^) 



1 



En 



2mNF{gA) 
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Pi ■P2 
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TO?, 



(4) 



where E2{Ei),Ep are, respectively, the initial (final) 
muon and the pion energies and Ai?i,. . — Ei,. — 
E^.. The functions F{gA) and G{gA) are given 

by: F{gA) = G{gA) - f^. One 

can easily check that Eq. ([4| is not factoriz- 
able into (production) x (propagation) x (detection) sub- 
processes due to the terms proportional to pi ■ P2 = 
E1E2 — \pi I \p2 1 cos a, where a is the angle between the di- 
rections of /x"*" particles. This is an important difference 
with respect to the case of neutrino-neutrino (AL = 0) 
oscillations where it was shown in Ref.[l2] that the S- 
matrix formalism reproduces the hypothesis of factoriza- 
tion of the probabilities. 

After integration over kinematical variables, it is pos- 
sible to write the rate of the complete process as 



'■Vi_^itE^ 



2E,, 



xF{M,(t>), (5) 



where F{M, (f) denotes the kinematical function 



F(M, 



2E, 



■{GFVudf\h\ M{gA~lf 



h - mNG{gA)h - 7:(™u - 



- 2mNF{gA) [mlh ~ (to^ + ml)h] ) . (6) 

The functions la for a — 1, .., 4 can be obtained from the 
following integral: 



Ia=lt^i^{E2- Ep)fJ{Ep + Em 



2E2 2E, 



■E 



N' 



'El 



-E2), 
(7) 
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with /i = 1, /2 = E2, fs = El, and /4 = {pi ■ P2) and tions. It started in 1982 when the BEBC bubble chamber 



_ (PrP2) 
J5 - E2-Ep ■ 

There are two interesting Umits to this process. At 
very short times, which means that the detection is very 
close to the production vertex (short-baseline neutrino 
experiment), one has, assuming that the E^. ~ E^, that 



El 



X F{M, i 



(8) 



where (m^^) is the effective Majorana mass for the muon 
neutrino. In the long time limit which corresponds to a 
long-baseline neutrino experiment, the oscillation terms 
cannot be neglected and this process depends on a new 
combination of phases, mixing angles and masses which 
could give us complementary information on the neutri- 
noless double beta decays or on any process that depends 
exclusively on the effective Majorana mass of the neutri- 
nos. Using this expression for the rate it is possible to 
get the CP asymmetry which will depend explicitly on 
Majorana phases. 



acp 



V - — r 

^i>] {U^iUy,iU*jU*j) m^^m^^ sin7 
T.^>3 Re {Uf,^Uf,^U;^U;^) m^^m^^ COS7 



(9) 
(10) 



where 7 — 



Here Ar 



is the difference 



2E^{GeV) ■ "'"'23 

in the squares of second and third eigenstate neutrino 
masses, Am^g = (2.43 ± 0.13) x lO'^ eV^, and L is 
the distance between production and detection vertices. 
Finally, Ei, is the energy of the neutrino beam. It is 
worth mentioning that the time evolution amplitude for 
the CP-conjugate process corresponds to the observa- 
tion oi fj,~ at the source and at the detector. There- 
fore, the associated nucleon weak vertex is given by 

(^Af'w)^ = UNiPNhi^igv +5A75)wjv'(Pjv')- In estimat- 
ing the CP asymmetry in Eq. pU]) , we have assumed that 

•^n'n — -^nn' ■ 

In the limit of ^13 — 0, the Majorana phases 
a I ^2 are the only sources of CP violation and hence 
lm{Uij,iUfj,iU*jU*j) cx sm{ai - aj). For i = 2 and j = 3 
one finds 



acp — tan [2(a2 - as)] sin 7. 



(11) 



Thus, in the case of long-baseline neutrino experiment 
like MINOS where the distance L is given by L = 735 
km and the energy E^ is typically around 2 — 3 GcV 
[3) [3 J one finds that sin 7 ~ 0{1). Thus, measuring 
CP asymmetry will be unavoidable indication for large 
CP violating Majorana phases. 



III. APPLICATION TO MINOS RESULTS ON 
NEUTRINO- ANTINEUTRINO OSCILLATIONS 

The last two decades have witnessed several experi- 
ments that investigate the neutrino-antineutrino transi- 



in the CERN SPS neutrino beam set a limit on — > 
and i>e thr oug h the search for v,, appearance. Re- 

cently, MINOS fl3| has measured the spectrum of 
events which are missing after travelling 735 km. It is 
these missing events which are the potential source of 
appearance. In their preliminary analysis, they were able 
to put a limit on the fraction of muon neutrinos transition 
to muon anti- neutrinos [ll|: 

P(i/^ - z>^) < 0.026 (90% C.I.). (12) 

Assuming CPT, this limit can be written as 



< 0.026 . 



(13) 



Using our expression for Pjy^-p,,, and the corresponding 
rate for neutrino oscillations Jl2], one gets 



< 0.001 



(14) 



In the limit of ultrarelativistic neutrinos, Ei,^ ~ ,£^(1 + 
m^. /2Ei,). and keeping the leading terms in the rriujEi, 
terms, we get 



< 0.001 X Et 



(15) 



To illustrate the usefulness of this relation, let us assume 
the case of two flavor neutrinos. In this case, one finds 

0.001X^2 > |(^^^)|2 

- 4 Re {UI2UII) m,,m,, sin" | 

- 2Im(C/22C/;|)m,,m,3sin7 . (16) 

From this equation, it is possible to get a bound on the 
effective muon-neutrino Majorana mass, only depending 
on the values of the Majorana phases as the oscillation 
terms cannot be neglected. If tt — 7/2 < 2(q2 — 0^3) < 
27r — 7/2 and using E^ k, 2 GeV, one gets the following 
conservative bound on 

l("^MM>l ;S64 MeV 

If —7/2 < 2(a2 — 03) < 0, the conservative bound on 
|(m^^)| is given by 

^ 109 MeV 

If < 2(a2 — 013) < TT — 7/2, it is not possible to get 



a conservative bound on | (r 



but the equation ([T 



could be used to bound Majorana parameters (masses 
and phases) which appear in |(m^^)|. 

So, within the hypothesis done on Majorana phases, 
the limits obtained improve by various order of magni- 
tude the actual limit on {man) coming from direct search 
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in t:~ ^'^ n'^ decay. Also, one could expect im- 

provement in experimental analysis in the next years and 
we could expect that with a specially designed neutrino 
experiment to measure these kind of processes, it should 
be possible to improve these limits by a few orders of 
magnitude in future experiments. 

In MINOS and in general in all long-baseline neutrino 
experiments, the oscillation terms are not negligible. So, 
it means that such analysis could not only give us in- 
formation on the neutrino effective Majorana mass but 
it could be used to determine parameters of the mixing 
matrices and bound on the absolute value of Majorana 
masses. Also, if in a long-baseline neutrino experiment, 
neutrino detectors are located at different distances from 
the source, it should be possible to get enough constraints 
on the mixing parameters and Majorana masses to fix 
them. 



IV. CONCLUSIONS 

If Majorana neutrinos do exist, | AL| ~ 2 processes like 
neutrino-antineutrino oscillations can occur. The pro- 
duction of leptons with same charges at the production 
and detection vertices of neutrinos will be a clear manifes- 
tation of these processes. In this paper we have used the 
S-matrix formalism of quantum field theory to describe 
these oscillations in the case of muon neutrinos produced 
in 7r+ decays which convert into muon antineutrinos that 
are detected via inverse beta decay on nucleons. 



One interesting result is that the time evolution prob- 
ability of the whole process is not factorizable into pro- 
duction, oscillation and detection probabilities, as is the 
case in neutrino oscillations We find that, for very 
short times of propagation of neutrinos, the observa- 
tion of /i^/i^ events would lead to a direct bound on 
the effective mass of muon Majorana neutrinos. In the 
case of long-baseline neutrino experiments, the CP rate 
asymmetry for production of fi'^ fi^ / fi~ fi~ events would 
lead to direct bounds on the difference of CP-violating 
Majorana phases. Finally, using the current bound on 
muon neutrino-antineutrino oscillations reported by the 
MINOS Collaboration we are able to set the bound 
(w^m) Si 64 MeV, which is several orders of magnitude 
below current bounds reported in the literature. 

As a consequence of these results, neutrino experi- 
ments aiming to measure neutrino-antineutrino oscilla- 
tions with different short- and long-baseline setups can 
be useful to get direct and complementary constraints on 
the masses and phases of Majorana neutrinos. 
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